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Abstract 

The  effect  of  solid  oxide  fuel  cell  (SOFC)  anode  thickness,  porosity,  pore  size,  and  pore  tortuosity  on  fuel  and  exhaust  gas  flow  is  calculated. 
Also  determined  is  the  concentration  of  these  gases  and  of  diluent  gases  as  a  function  of  position  across  the  anode.  The  calculation  is  based  on 
the  dusty-gas  model  which  includes  a  Knudsen  (molecule-wall)  collision  term  in  the  Stefan-Maxwell  equation  which  is  based  on  unlike-molecule 
collisions.  Commonly  made  approximations  are  avoided  in  order  to  obtain  more  exact  results.  One  such  approximation  is  the  assumption  of  uniform 
total  gas  pressure  across  the  anode.  Another  such  approximation  is  the  assumption  of  zero  fuel  gas  concentration  at  the  anode-electrolyte  interface 
under  the  anode  saturation  condition  for  which  the  SOFC  output  voltage  goes  to  zero.  Elimination  of  this  approximation  requires  use  of  a  model  we 
developed  (published  elsewhere)  for  terminal  voltage  V  as  a  function  of  electrolyte  current  density  i.  Key  formulae  from  this  model  are  presented. 
The  formulae  developed  herein  for  gas  flow  and  tortuosity  are  applied  to  the  results  of  a  series  of  careful  experiments  performed  by  another  group, 
who  used  binary  and  ternary  gas  mixtures  on  the  anode  side  of  an  SOFC.  Our  values  for  tortuosity  are  in  a  physically  reasonable  low  range,  from  1 .7 
to  3.3.  They  are  in  fair  agreement  with  those  obtained  by  the  other  group,  once  a  difference  in  nomenclature  is  taken  into  account.  This  difference 
consists  in  their  definition  of  tortuosity  being  what  some  call  tortuosity  factor,  which  is  the  square  of  what  we  and  some  others  call  tortuosity.  The 
results  emphasize  the  need  for  careful  design  of  anode  pore  structures,  especially  in  anode- supported  SOFCs  which  require  thicker  anodes. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  provide  efficient  generation  of  electricity  using  a  variety  of  fuels.  They  produce  only  steam  exhaust 
if  hydrogen  is  the  fuel.  The  same  device  can  act  as  a  steam  electrolyzer  to  produce  hydrogen  and  oxygen  from  steam  if  a  dc  voltage 
is  applied. 

Present  SOFCs  operate  in  the  800-1000  °C  (1073-1273  K)  range,  but  strong  efforts  are  in  progress  to  extend  this  range  downward 
to  600  °C  (873  K)  or  lower  to  minimize  interconnect  cost  by  allowing  use  of  stainless  steel,  and  to  reduce  problems  associated  with 
thermal  expansion,  atomic  migration,  and  corrosion.  The  main  obstacle  for  this  temperature  reduction  is  the  rapid  decrease  in  solid 
electrolyte  oxygen  ion  conductivity  with  decreasing  temperature.  Accordingly,  the  electrolyte  is  made  as  thin  as  possible,  of  order 
10  fxm.  This  requires  that  an  electrode  layer,  usually  the  anode,  must  be  made  thick  enough  to  mechanically  support  the  cell. 

With  a  thicker  anode,  one  must  calculate  carefully  the  fuel  and  exhaust  gas  concentrations  at  the  anode-electrolyte  interface. 
These  concentrations  (or  corresponding  partial  pressures)  strongly  affect  V(i),  the  terminal  voltage  V  as  a  function  of  electrolyte 
current  density  i.  The  anode  characteristics  (assuming  a  much  thinner  cathode)  also  under  some  conditions  determine  the  maximum 
attainable  i.  This  is  the  current  density  when  the  fuel  flow  rate  is  high  enough  to  reduce  the  fuel  gas  partial  pressure  at  the  interface 
to  near  zero,  and  the  SOFC  output  voltage  to  exactly  zero.  This  fuel  pressure  gradient  in  the  anode  is  increased  by  the  inflowing  fuel 
molecule  collisions  with  the  outflowing  exhaust  molecules. 
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Table  1 

Tortuosities  calculated  for  H2/H2O  and  CO/CO2  binary  systems 


Fuel 

p.p.  (atm) 

/as  (A  cm *  2) 

r 

r  (ftfas=0) 

(t,v)1/2 

Tjv 

h2 

0.2 

1.3 

1.945 

2.179 

2.211 

4.89 

h2 

0.34 

2.1 

2.024 

2.243 

2.294 

5.26 

h2 

0.5 

3.25 

1.913 

2.136 

2.211 

4.89 

CO 

0.18 

0.21 

2.730 

2.871 

2.876 

8.27 

CO 

0.23 

0.31 

2.502 

2.656 

2.661 

7.08 

CO 

0.32 

0.49 

2.311 

2.475 

2.484 

6.17 

CO 

0.44 

0.73 

2.194 

2.366 

2.377 

5.65 

Fuel  partial  pressures  (p.p.),  anode  saturation  current  densities  (/ 

as),  and  tjv  are  from  JV  Table  2  [1]. 

Table  2 

Tortuosities  calculated  for  various  ternary  system  diluent  gases  and  partial  pressures  (p.p.  in  atm),  and  measured  anode  saturation  current  densities  (/as,  A  cm 

-2) 

Diluent 

p.p. 

las 

Hill 

”^quad 

T exact 

(Dv)1/2 

Trv 

He 

0.78 

1.5 

1.7159 

1.7025 

1.7025 

2.3707 

5.62 

He 

0.65 

2.4 

1.7295 

1.7069 

1.7068 

2.3130 

5.35 

He 

0.53 

3.25 

1.7371 

1.7052 

1.7051 

2.2583 

5.1 

He 

0.42 

4 

1.7634 

1.7215 

1.7214 

2.2271 

4.96 

n2 

0.8 

0.89 

2.8705 

2.8020 

2.8029 

2.7495 

7.56 

n2 

0.67 

1.56 

2.7149 

2.6225 

2.6240 

2.6627 

7.09 

n2 

0.57 

2.3 

2.4401 

2.3469 

2.3484 

2.4597 

6.05 

n2 

0.5 

2.65 

2.4467 

2.3433 

2.3450 

2.4860 

6.18 

co2 

0.81 

0.7 

3.4465 

3.3183 

3.3225 

3.0000 

9 

co2 

0.68 

1.43 

2.9452 

2.8038 

2.8094 

2.6833 

7.2 

co2 

0.6 

1.88 

2.8036 

2.6560 

2.6622 

2.6153 

6.84 

co2 

0.5 

2.4 

2.7221 

2.5663 

2.5730 

2.6019 

6.77 

The  first  three  and  last  columns  are  from  JV  [1]. 


A  key  parameter  in  determining  the  pressure  gradient  is  the  tortuosity  r,  which  is  the  ratio  of  the  typical  diffusion  path  length 
to  the  electrode  thickness.  Measurement  of  its  value  is  needed  for  determining  the  quality  of  the  anode  pore  configuration.  Also, 
t  is  a  necessary  parameter  in  any  model  for  V(i),  and  conversely  we  will  show  how  our  model  for  V(i)  provides  a  more  accurate 
determination  of  t. 

The  paper  begins  with  a  description  of  the  Jiang-Virkar  [1]  experiment  to  which  our  analysis  is  applied.  This  is  followed  by  a 
discussion  of  the  Stefan-Maxwell  equation  including  both  Knudsen  (molecule  collisions  with  pore  walls)  and  binary  collision  effects. 
Then  the  theory  and  results  for  binary  and  ternary  anode  gas  mixtures  are  presented,  followed  by  conclusions  and  recommendations. 


2.  Description  of  Jiang-Virkar  experiment 


Jiang  and  Virkar  [1]  performed  a  series  of  careful  experiments  on  flow  of  various  gas  mixtures  in  the  anode  of  an  SOFC.  This 
paper  will  be  referred  to  as  JV.  JV  studied  binary  system  anode-side  gas  mixtures  H2/H2O  and  CO/CO2,  and  ternary  systems 
with  incoming  H2  fuel  diluted  with  He,  N2,  or  CO2.  In  all  cases,  the  total  metered  inflow  rate  of  the  two  gases  was  held  at 
140  mL  min-1 ,  corresponding  to  6.286  x  10 19  molec  s-1  as  explained  below.  The  “ternary”  systems  are  called  ternary  because  H2O 
(steam)  “exhaust”  is  produced  at  the  anode-electrolyte  interface  and  mixes  with  the  hydrogen  and  diluent  gas.  We  compare  their 
calculations  of  tortuosity  r  with  values  found  from  our  analysis.  We  use  the  definition  of  tortuosity  as  the  average  actual  path  length 
traversed  through  the  anode  by  a  fuel  or  exhaust  molecule,  divided  by  the  anode  thickness.  There  are  other  definitions  in  the  literature; 
in  fact,  the  term  “tortuosity”  has  a  “tortuous”  history  [2].  In  particular,  JV  and  some  others  define  tortuosity  as  what  we  and  some 
others  [2]  call  “tortuosity  factor.”  The  low  range  of  r  values  we  obtain,  from  1.7  to  3.3,  is  in  accord  with  the  2.5-3  range  found  by 
Williford  et  al.  [3]  who  criticized  reports  of  “anomalously  high”  tortuosities  in  the  6-17  range. 

In  the  JV  experiments,  the  cathode- side  gas  was  air.  Both  sides  must  have  been  at  1  atm  sea  level  in  order  to  reproduce  their 
results,  even  though  their  laboratory  is  at  the  University  of  Utah,  well  above  sea  level.  The  partial  pressures  in  atmospheres  for  the 
various  metered  input  gases  are  given  in  our  Tables  1  and  2  and  Figs.  2  and  3.  The  temperature  T  was  800  °C  =  1073  K  for  all  their 
experiments.  Their  SOFC  consisted  of  five  layers,  as  shown  in  Fig.  1,  and  constructed  of  the  materials  listed  below. 


Anode 

Anode  interlayer 
Electrolyte 
Cathode  interlayer 
Cathode 


Ni  + YSZ  (yttria- stabilized  zirconia) 

Ni  +  YSZ 

YSZ  +  SDC  (samaria-doped  ceria)  bilayer 

LSC  (Sr-doped  LaCoC>3)  +  SDC 

LSC 
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Cathode:  LSC 

Cathode  interlayer:  I.SC  (Sr-doped  I-aCoO,)  +  SIX' 


Anode  interlayer:  Nl  +  YS7. 

Anode:  Ni  +  YSZ  (yttria -stabilized  zirconia) 


Fig.  1.  Drawing  of  SOFC  employed  by  Jiang  and  Virkar  [1], 


The  anode  thickness  w  =  0.11  cm,  and  its  porosity  0  =  0.54.  Their  estimated  pore  radius  from  SEM  micrographs  was 
-0.5  [Jim. 


3.  Dusty-gas  model  for  gas  flow  in  pores 


The  dusty-gas  model  consists  of  the  Stefan-Maxwell  equation  including  the  Knudsen  term.  In  molar  units  it  is  [1] 
Ni  ^ XjNi-XiNj  (1  /RT)d(PXi) 

—  = - S — • 


where  Ni  and  Nj  are  molar  fluxes  of  components  i  and  j  (mol  cm-2  s_1),  respectively,  Dk/  and  D/y  are  the  Knudsen  diffusion 
coefficient  for  component  i  and  the  binary  diffusion  coefficient  for  components  i  and  j,  respectively,  Xi  and  Xj  are  the  fractional 
molar  concentrations  of  components  i  and  /,  P  is  the  total  pressure,  R  is  the  gas  constant,  and  T  is  the  absolute  temperature.  We 
emphasize  that  x  is  the  coordinate  along  a  typical  diffusion  path,  not  across  the  anode,  so  the  total  diffusion  path  length  is  zw ,  the 
product  of  the  tortuosity  and  the  anode  thickness.  Note  that  this  analysis  ignores  laminar  flow  (permeation  flux)  effects  [4]  which 
are  considered  unimportant  in  small  pores. 

Note  further  that  Eq.  (1)  is  modified  from  Eq.  (9)  of  the  JV  paper  in  that  P  is  outside  the  derivative  in  their  paper.  That  is  correct 
only  in  the  approximation  that  P  is  independent  of  position.  The  total  pressure  of  the  incoming  gases  and  the  gases  in  the  plenums 
is  one  atmosphere,  denoted  as  P\ ,  but  the  total  pressure  P  in  the  anode  and  cathode  pores  will  be  a  function  of  position  in  the  pore. 
We  will  show  that  P  varies  across  the  anode  by  amounts  of  order  10%. 

The  Knudsen  diffusion  coefficient  only  considers  collisions  with  the  wall,  and  from  the  kinetic  theory  of  gases  is 


Dkj  = 


2  /8 RT\l/2_ 
3 \nM~i)  ' 


(2) 


where  r  is  the  mean  pore  radius  (0.5  |mm)  and  Mi  is  the  molar  mass.  At  1073  K  the  £>kj  values  for  this  pore  radius  are  11.3  and 
3.00  cm2  s_1  for  H2  and  CO,  respectively.  To  convert  Eq.  (2)  to  molecular  units,  we  replace  R/Mi  by  k/m;,  where  k  is  Boltzmann’s 
constant,  1.38  x  10-23  J  K-1  and  m;  is  molecular  mass.  We  convert  to  molecular  units  to  make  some  equations  simpler. 

The  binary  diffusion  coefficients  Dy  for  various  gas-molecule  pairs  at  three  temperatures  are  tabulated  in  JV  Table  1.  The 
coefficients  for  the  two  2-species  cases,  H2/H2O  and  CO/CO2,  have  Dy\  values  of  7.704  and  1.408  cm2  s-1,  respectively  at  1073  K. 
The  subscript  1  indicates  that  the  values  are  for  1  atm  total  pressure. 

To  convert  Eq.  (1)  to  molecular  units,  several  steps  are  necessary.  The  Xi  are  dimensionless  mole  or  molecular  fractions  of  the 
total  gas  concentration.  The  Ni  are  converted  to  Ji  molecules  per  unit  area  and  time  by  multiplying  by  Avogadro’s  number  Na-  The 
gas  constant  R  has  moles  in  the  denominator,  so  putting  N\  in  the  numerator  changes  R  to  k.  But  P  =  nkT ,  where  n  is  the  total  number 
of  gas  molecules  per  unit  volume,  so  the  kT  factors  cancel.  Since  nXi-nu  the  right  side  of  Eq.  (1)  is  simply  —drii/dx.  Dy  has  the 
total  concentration  n  in  its  denominator,  but  since  Dy  itself  is  in  the  denominator  in  Eq.  (1),  n  appears  in  the  numerator  in  that  term, 
and  n  multiplied  by  Xi  or  Xj  changes  them  to  nt  or  nj.  Thus,  Eq.  (1)  converts  to 

Ji  (Jjjtj  -  JjUi)kT  dm 

DKi  DijiPi  dx'  U 

Here,  Dy\  is  the  binary  diffusion  coefficient  at  P\  =  1  atm,  and  the  nkT- type  terms  in  the  numerator  divided  by  P\  convert  Dy\  to 
Dy  at  the  actual  total  pressure  at  any  position  x  along  the  anode  pore. 

The  variable  v  is  distance  along  the  pore,  so  the  total  length  L  of  the  pore  is  L  =  zw,  where  z  is  the  anode  pore  tortuosity  and  w 
is  the  0.1 1  cm  thickness  of  the  anode. 
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4.  Calculation  for  binary  gas  input 


For  the  two  binary  gas  systems  H2/H2O  and  CO/CO2,  if  subscript  f  is  for  the  fuel  gas  and  e  is  for  the  exhaust  gas,  /e  =  —Jf 
in  Eq.  (3).  Here,  Jf  =  iz/4>q,  where  i  is  the  current  density  (A  cm-2  =  Ccm-2s-1)  in  the  solid  electrolyte  and  q  is  the  charge 
(3.2  x  10“ 19  C  molecule-1)  carried  per  gas  molecule  annihilated/created  in  the  reaction  at  the  anode/electrolyte  interface. 

The  r/0  factor  is  an  enhancement  factor  by  which  the  flow  density  Jt  is  enhanced  compared  to  its  value  if  the  anode  were 
completely  porous  (0  =  r  =  1).  To  derive  this  result,  note  that  this  enhancement  factor  for  an  anode  of  area  S  with  Apa  pores  each 
of  cross-sectional  area  Ax  is  S/NV3iAx.  The  porosity  0  is  the  total  pore  volume  NvaAxzw  divided  by  the  anode  volume  Sw ,  so 
0  =  NVSiAxz/S.  By  rearranging  this  equation  we  find  that  the  enhancement  factor  S/NvslAx  is  r/0. 

Eq.  (3)  for  the  fuel  gas  can  now  be  written,  after  setting  nQ  +  rif  =  n ,  as 


dHl  =  _(iL'\  (2-4-  =  _r,r  _ 


dx 


<Pq)  V^Kf  +  D&nPi ) 


C  3T  —  C4HT , 


(4) 


where  C3  and  C4  are  positive  constants,  and  subscript  f  refers  to  the  fuel  gas  (H2  or  CO)  in  the  binary  gas  system.  The  corresponding 
equation  for  the  exhaust  gas,  after  defining  (mjrrif)112  =  2h  +  1,  is 


3  ne 
dx 


=  (2  h  +  1)C3T  +  C4UT. 


Adding  Eqs.  (4)  and  (5)  yields  the  equation  for  n  and  its  solution, 
3  n 

—  =  2hc?>x,  n  =  n\  +  2hc?>zx. 
dx 


(5) 


(6) 


From  the  ideal  gas  law,  P  =  nkT ,  we  see  that  the  total  pressure  P  increases  linearly  with  x,  in  going  through  the  anode  from  the 
plenum  to  the  solid  electrolyte. 

Now  that  n(x)  is  known,  we  can  solve  Eq.  (4)  for  rif. 


—  =  —c^z  —  C4(n\  +  2hc^zx)z, 
dx 


rif  =  ftfp  —  (C3  +  C4n\)zx  —  hc?,C4Z2x2, 


(7) 


where  rifv  is  the  concentration  of  fuel  gas  in  the  plenum,  which  will  be  found  below.  We  see  that  rif ,  and  thus  nQ  =  n  —  rif ,  vary  both 
linearly  and  quadratically  with  x. 

If  the  cell  current  density  is  limited  by  anode  saturation,  designated  as  4s  by  JV,  then  i  becomes  4s  in  Eq.  (4)  and  in  the  constants 
C3  and  C4,  which  are  designated  C3as  and  C4as  for  the  anode  saturation  condition.  In  previous  works  [1,5]  the  approximation  was 
made  that  in  this  limit  the  fuel  concentration  rif  reaches  zero  at  the  anode-electrolyte  interface  where  x  =  zw.  This  approximation 
obviously  cannot  be  exactly  correct,  because  4s  cannot  be  positive  for  zero  fuel  concentration,  whereas  the  measured  4s  is  positive. 

To  avoid  making  this  unphysical  approximation,  we  use  results  from  our  recent  analysis  of  SOFC  voltage  V  as  a  function  of 
current  density  i  [8].  This  relation  is 


m  = 


U  kT 

- In 

q  q 


(b  +  iXd  +  iff 

(a  —  i)(c  —  z')1/2 


Vohm. 


(8) 


Here,  U is  the  enthalpy  for  the  gas-state  reaction  H2+I/2O2  ->  H2O  or  CO  +  I/2O2  ->  CO2  at  800  °C  and  1  atm  for  all  constituents, 
q  is  the  reaction  charge  transfer  of  2e,  k  is  Boltzmann’s  constant,  T  is  absolute  temperature,  and  Vohm  is  the  ohmic  polarization. 
Unlike  other  V(0  expressions  in  the  literature  (see  [8]),  Eq.  (8)  puts  all  four  reaction  attempt  current  densities  a ,  b,  c ,  d  on  the  same 
footing.  The  forward  and  reverse  anode  attempt  current  densities  are  a  and  b ,  while  the  forward  and  reverse  cathode  attempt  current 
densities  are  c  and  d.  These  reaction  attempt  current  densities  are  multiplied  by  reaction  success  probabilities  which  are  functions  of 
k ,  T,  q ,  /,  and  anode  and  cathode  contributions  to  U ,  to  find  anode  and  cathode  contributions  to  V  in  Eq.  (8).  Eq.  (15)  shows  the  form 
of  a ,  and  the  other  incoming  gas  reaction  attempt  current  densities  b  and  c  have  similar  forms.  Other  authors  (see  [8])  have  factors 
in  the  In  argument  corresponding  more  or  less  to  ours  that  involve  a ,  b ,  and  c,  though  these  factors  involve  partial  pressures  whereas 
ours  involve  concentrations  and  molecular  velocities.  No  other  authors  to  our  knowledge  explicitly  show  a  factor  corresponding  to 
our  (d+  i)m  factor,  where  d  is  the  attempt  current  density  corresponding  to  oxygen  ions  from  the  electrolyte  combining  to  form  O2 
gas.  This  factor  is  much  larger  than  the  other  three  factors  in  the  In  argument,  and  reduces  the  open-circuit  emf  considerably  below 
the  value  of  1.214  V  for  Ulq.  Other  authors  instead  implicitly  incorporate  the  —(kT/q)  \n[(d  +  /)1/2]  term  as  an  entropy  term  as  part 
of  the  reaction  Gibbs  energy  G,  and  present  the  known  G  in  Glq  in  place  of  our  Ulq  as  the  first  term  in  V(i).  Fortunately,  our  d 
expression  [8]  based  on  a  molecular  model  for  that  reaction  attempt  rate  gives  approximately  the  correct  open-circuit  emf  and  thus 
provides  a  check  on  that  reaction  model. 

The  terminal  voltage  V  is  often  expressed  as  a  sum  of  terms 


E  —  Vq  Vact  VCOnc  E0hm> 


(9) 
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and  each  term  can  be  found  separately  as  shown  below.  The  open-circuit  emf  Vo,  from  Eq.  (8)  with  i  =  Vohm  =  0,  and  with  a ,  b ,  c,  d 
set  at  their  open-circuit  values  ao,  bo,  cq,  do,  is 


V0  = 


U  kT  bod!/2 

—  —  —  In - - — 

A  n  1/2  * 

<1  a0c0 


(10) 


The  activation  polarization  from  Eq.  (9)  is  given  by  Vact  =  Vo  —  V—  VCOnc  —  V0hm-  It  signifies  the  voltage  drop  required  to  provide 
a  given  i  for  fixed  gas  concentrations  at  the  electrode  interfaces  with  the  electrolyte.  If  we  set  a  =  ao ,  etc.  in  Eq.  (8)  we  eliminate 
Vconc  because  the  gas  concentrations  at  the  electrode  interfaces  with  the  electrolyte  depend  only  on  a ,  b ,  and  c.  If  we  also  eliminate 
Vohm  in  Eq*  (8),  then  Eq.  (8)  becomes 


Vo  -  Vact 


U  kT  (bo  +  Dido  +  i)1/2 

- In - Yn  ’ 

q  q  (at)  -  i)(c0  -  i) 1 


(11) 


so  Vact  is  found  by  subtracting  the  Eq.  (11)  from  the  Eq.  (10)  expression, 

kT  aoc\/2(bo  +  i)(do  +  i)l/2 
Vact  =  —  In  — ^ 

<7  b0df  l(ao  -  i)(co  -  i)'/2 


Inserting  expressions  from  Eqs.  (8)  and  (12)  into  Eq.  (9),  we  obtain 

_kT  (ao-  i)(b  +  Q(cq  -  i)l/2(d  +  z)1/2 
C°nC  q  (a  -  i)(bo  +  i)(c  -  /)1/2(do  +  /)1/2 


The  meaning  of  Vconc  is  the  additional  voltage  drop  required  because  increasing  i  causes  changes  in  gas  concentrations  at  the 
electrode  interfaces  with  the  electrolyte.  An  important  fact,  generally  not  mentioned,  is  that  there  are  two  contributions  to  VCOnc-  The 
one  that  usually  comes  to  mind  is  the  change  with  i  of  gas  concentrations  as  one  goes  across  the  electrodes.  The  other  contribution 
is  the  dependence  on  i  of  gas  partial  pressures  and  concentrations  in  the  plenums  which  occurs  even  if  there  is  no  concentration 
change  across  the  electrodes  and  the  feed  gas  flow  rates  remain  fixed.  Both  contributions  are  discussed  below.  Both  Vact  and  VCOnc 
take  on  the  sign  of  i.  The  sum  of  Vact  and  Vconc  is 


Vact  +  Vconc 


kT  aocl0/2(b  +  i)(d  +  /)1/2 
q  bod^2(a  —  i)(c  —  z)1/2 


(14) 


The  saturation  current  density  zas  is  defined  as  the  current  density  at  which  V  drops  to  zero,  but  from  both  experiment  and  our 
V(z)  model  [8]  the  voltage  drops  extremely  rapidly  toward  —  oo  as  i  goes  through  zas.  Accordingly,  to  a  very  good  approximation,  zas 
is  the  current  density  at  which  (an  applied)  V  approaches  —  oo.  From  Eq.  (8),  this  occurs  for  zas  =  a. 

Our  model  [8]  gives  a  relation  between  a  and  zzf,  allowing  us  to  find  the  nf  value  zzfas  attained  at  the  anode  saturation  condition 
where  a  =  zas.  This  relation  is 


1 


a  =  -n  f 
2 


1/2 

qi\  -  Ua)/a, 


(15) 


where  nif  is  the  fuel  molecule  mass,  (j kT/rrif )1/2  is  the  fuel  molecule  average  velocity  component  (through  the  anode  pore)  directed 
toward  the  triple  phase  boundary  (TPB),  1  —  va  is  the  probability  that  an  O2-  site  at  the  anode  TPB  will  be  occupied  (not  vacant)  so 
the  fuel  oxidation  reaction  can  occur,  and /a  is  the  probability  that  the  fuel  molecule  strikes  the  electrolyte  at  an  anode  TPB.  The  Eq. 
(15)  relation  is  model-dependent,  but  it  seems  that  in  any  model  there  will  be  a  linear  or  nearly  linear  relation  between  a  and  rif. 

Going  back  now  to  Eq.  (7),  we  can  no  longer  set  rif(x  —  L  —  zw)  =  0,  but  instead  we  designate  it  as  zzfas.  We  then  obtain  a 
quadratic  equation  for  t2,  valid  for  the  anode  saturation  condition 

^C3aSC4asW2r4  +  (C3as  +  C4as«l)WT2  -  Wfpas  +  «fas  =  At4  +  Bz2  +  C  =  0.  (16) 


All  constants  in  Eq.  (16)  except  C  are  positive,  so  in  the  solution 
,  —B  ±  (B2  -  4 AC)1/2 

T  =  - , 

2A 


(17) 


we  must  choose  the  +  sign  so  that  r2  will  be  positive. 

All  constants  in  Eq.  (16)  have  been  defined  except  for  the  boundary  conditions  n\  and  zzfpas  for  the  total  and  fuel  gas  concentrations, 
respectively  in  the  plenum  at  the  outer  anode  surface,  and  zzfas  for  the  fuel  gas  concentration  at  the  anode-electrolyte  interface.  For 
the  JV  experiments,  n\  =P\/kT  where  Pi  =  1  atm=  1.015  x  105  N m-2  and  T  =  1073  K,  so  n\  =  6.855  x  1018  cm-3. 
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To  find  ftfpas  we  consider  the  total  flow  in  and  out  of  the  anode  plenum.  From  the  discussion  in  JV  p.  A943  it  appears  that  for  all 
the  binary  and  ternary  systems  displayed  in  JV  Table  2,  the  metered  gas  inflow  rate  of  140  mL  min-1  is  for  two  gases,  so  that  in  the 
ternary  systems  the  third  gas  is  H2O  generated  at  the  anode-electrolyte  interface.  This  total  metered  inflow  rate  of  140  mL  min-1  at 
a  temperature  of  T-m  =  273  K  and  at  1  sea  level  atmosphere  is,  from  the  ideal  gas  law, 


Pi 


dVn 
d  t 


kT'm 


dNw 


d  t 


so 


dNin 

dt 


=  ji  =  6.286  x  1019molecs  1 . 


(18) 


JV  did  not  specify  the  temperature  at  the  metered  inflow,  but  Tm  =  273  K  fits  their  results  best,  based  on  solving  for  mj  from 
Eq.  (32)  of  JV  using  £>H2,eff>  4is,  and  pm  from  JV  Table  2  for  the  three  H2/H2O  runs.  The  metered  inflow  of  each  gas  constituent 
is  found  by  multiplying  the  total  metered  inflow  rate  by  the  partial  pressure  fraction  given  or  implied  in  JV  Table  2.  For  instance, 
the  first  entry  is  pf  =  0.2  for  H2,  so  pQ  =  0. 8  for  H2O.  Though  these  Table  2  entries  are  in  atmospheres,  in  this  analysis  we  define 
these ps  as  dimensionless  partial  pressure  fractions.  Then  the  respective  metered  inflow  rates  are:  jfm=Pfji  =  1-257  x  1019  s_1  and 
jem  =  Peji  =  5.029  x  1019  s~l .  We  assume  here  that  the p  values  in  JV  Table  2  are  really  p°  values,  where p°  is  the  partial  pressure  of 
the  metered  incoming  fuel  or  diluent  as  defined  below  Eq.  (26)  of  JV. 

To  find  the  net  inflow  of  each  gas  into  the  plenum,  we  must  also  consider  the  fuel  gas  outflow  into  the  anode,  and  the  exhaust 
gas  inflow  from  the  anode.  These  flows  are  equal  and  opposite,  and  are  given  by  ye  =  —jf  =  iS/q ,  where  S  is  the  active  area  of  1.1  cm2 
and  q  =  3.2  x  10-19  C.  Using  again  the  first  entry  in  JV  Table  2,  zas  =  1.3  A  cm-2,  we  find  H2  outflow y'fas  =  —4.469  x  1018  molec  s-1 
and  H2O  inflow  7eas  -+4.469  x  1018  molec  s-1.  Comparison  of  the  H2  outflow  with  the  metered  H2  inflow  gives  a  fuel  utilization 
— j'fas/jfm  of  35.5%  for  this  example. 

From  the  numbers  in  the  previous  two  paragraphs,  the  net  H2  inflow  into  the  plenum  is  jfnQt  =  7fm  +  7fas  =  8.104  x  1018  molec  s-1 
and  the  net  H2O  inflow  is yenet  =  jem  +  7e as  =  5.476  x  1019  molec  s_1 ,  for  a  total  net  inflow  to  the  plenum  of  j\  =  6.286  x  1019  s-1.  The 
outflow  from  the  plenum  to  the  atmosphere  must  contain  the  corresponding  fractions,  — /fnet//l  =  0.1289  for  H2  and y'enet/ji  =  0.871 1 
for  H2O,  on  a  molecular  basis.  These  values  are  considerably  different  from  the  corresponding  metered  inflow  fractions,  pf  =  0.2 
and  pQ  -  0.8.  To  summarize  the  analytic  expression  for  jfnet  and yenet> 


iS  .  iS 

Jfnet  =  PfJ  1 - ;  Jenet  =  (1  -  Pf)jl  +  —  - 

q  q 


(19) 


One  must  consider  whether  the  mole  fractions  in  the  plenum  are  the  same  as  the  mole  fractions  in  the  outflow.  For  openings  small 
compared  to  the  mean  free  path,  this  is  not  the  case,  because  the  outflow  rate  per  molecule  in  the  plenum  is  inversely  proportional 
to  the  square  root  of  the  molecular  mass  [6].  However,  it  seems  probable  that  the  outflow  arrangement  provided  by  JV  result  in  the 
mole  fractions  in  the  plenum  being  essentially  the  same  as  in  the  outflow.  We  assume  that  the  gases  in  the  plenum  are  well  mixed  and 
have  the  same  mole  fractions  everywhere.  Accordingly,  in  the  plenum  the  fuel  and  exhaust  concentrations  under  anode  saturation 
conditions  are 


^fpas 


1  -  Pf  + 


»1, 


(20) 


where  pr  and  zas  are  values  taken  from  JV  Table  2,  and  i\  =  /)  q/S  =  1 8.29  A  cm-2. 

The  remaining  parameter  to  determine  is  /ifas,  the  fuel  gas  concentration  at  the  anode-electrolyte  interface  under  the  anode 
saturation  condition.  From  Eq.  ( 1 5),  rif  remains  proportional  to  a  as  i  changes,  because  we  assume  the  other  parameters  are  independent 
of  i.  For  open-circuit  (i  =  0)  conditions,  rif  at  the  interface  equals  rifo  at  the  plenum,  which  from  Eq.  (20)  will  be  pfU\  for  open- 
circuit.  This  known  value  for  open-circuit  rif  inserted  into  Eq.  (15)  provides  a  known  value  ao  =pfao\  for  the  open-circuit  value 
of  a ,  where  aoi  is  the  open-circuit  value  for  pf=l.  For  the  anode  saturation  condition,  a  =  /as,  and  because  rif  is  proportional  to 
a ,  nfSLS  =  rifQaSLS/ao=pfniiSLS/pfaoi  =  n\idiJaf)\.  Inserting  this  rifas  value  and  the  above  ftfpas  value  into  Eq.  (16)  yields  the  following 
equation  for  t2  in  terms  of  known  parameters 


^3asc4asu;2T4  +  (c3as  +  C4asni)wT2  - 


n\  +  n\ -^-  =  At4  +  Bz2  +  C  =  0. 
<201 


In  our  previous  analysis  [5]  the  ^i/as/^oi  term  was  missing. 


(21) 


5.  Results  for  binary  gas  input 


We  are  now  in  a  position  to  use  Eqs.  (21)  and  (17)  to  find  z 


f  -B  +  (B2  -4AC) 


1/2 


1/2 


Z  = 


2  A 


(22) 
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Fig.  2.  Plots  of  gas  concentrations  vs.  position  in  anode  for  some  binary  gas  inputs  under  anode  limiting  current  conditions. 


We  begin  by  summarizing  the  constants  that  appear  directly  or  indirectly  in  Eq.  (16),  and  that  do  not  change  from  run  to  run. 
They  are  provided  here  to  aid  the  reader  who  may  want  to  check  the  numerical  results  of  this  study  that  appear  in  Table  1  and  Fig.  2. 
The  meanings  of  the  symbols  are  explained  above  where  the  symbols  first  appear. 


0  =  0.54,  4  =  3.2  x  1(T19C,  k  =  1.38  x  10“ZJJK“\  P=  1073K,  Pi  =  1.015  x  103 Nm“z, 

u>  =  0.11cm,  n\  =  6.855  x  1018 cm-3,  i\  =  18.29Acm-2,  am  =  48.4 Acm-2, 

Z?Kf  =  11.3cm2s_1(H2);  3.00  cm2  s-1  (CO), 

£>efl  =  7.704cm2  s-1  (H2,H20)=  1.408  cm2  s-1(CO,  C02), 

—  =  — - —  =  5.121  x  1017sC-1  cm-2(H2)  =  1.929  x  1018sC-1  cm-2(CO), 

('as  (pqDjtf 

^  =  — — - =  0. 1096  cm  s  C-1  (H2)  =  0.5994cmsC-1(CO), 

*as  0*7  f^ef  1 P 1 

h  =  (We/m f)  ~  1  =  1(H2),  =  0.1268(CO), 


-23  T  TZ-1 


(23) 


/  Zas  ias  \  ^ 

ftfpas  =  -C  =  [Pi  -  - - n  1,  A  =  hc3asC4asW  ,  B  =  (c3as  +  C4aSni)w. 

V  n  aoi  J 

We  carry  out  a  sample  calculation  for  the  first  entry  in  JV  Table  2,  for  the  H2/H2O  binary  system,  for  which  /?f  =  0.2  and 
zas  =  1.3  Acm-2.  Various  parameters  used  in  this  calculation  and  whose  meanings  are  explained  above  are 

qqi  =48. 4 Acm-2,  ao  =  pfaoi  =9. 68 Acm-2  from  our  V(i)  analysis  [8] 

c3as  =  6.657  x  1017  cm-4,  c4as  =  0.1425  cm-1,  A  =  1.148  x  1015cm-3,  B  =  1.807  x  1017cm-3, 

C  =  -6.996  x  1017  cm-3  (Ref.  [5]  value  -  8.836  x  1017  cm-3),  (24) 

(. B 2  -  4AC)1/2  =  1.894  x  1017  cm-3  (Ref.  [5]  value  1.916  x  1017  cm-3), 
r  =  1.944  (Ref.  [5]  value  2.179),  JVr1/2  =  2.211,  JVr  =  4.89. 


These  results  and  those  for  the  other  two  H2/H2O  data  sets  from  JV  appear  in  Table  1. 

It  is  highly  encouraging  that  our  r  values  agree  well  with  each  other,  and  especially  that  their  small  magnitudes  seem  very 
reasonable  for  the  high-porosity  (0  =  0.54)  anode  that  JV  investigated.  They  agree  fairly  well  also  with  the  JV  results  in  the  (rjy)1/2 
column,  where  it  is  taken  into  account  that  what  JV  call  tortuosity  in  the  rjy  column  is  the  tortuosity  factor  [2]  whose  square  root  is 
what  we  and  others  call  tortuosity.  To  calculate  the  r  values  for  the  four  data  sets  from  JV  for  the  CO/CO2  binary  system,  the  above 
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analysis  is  unchanged  except  that  subscripts  f  and  e  refer  to  CO  and  CO2,  respectively.  The  tortuosity  results  are  given  in  Table  1 
below. 

We  see  in  Table  1  that  our  previous  [5]  analysis  results  [r(ftfas  =  0)]  which  range  from  about  2.1  to  2.9  agree  better  than  our 
present  results  [t]  which  range  from  about  1.9  to  2.7,  with  the  Jiang- Virkar  [1]  results  [(rjy)1/2]  which  range  from  about  2.2  to 
2.9.  This  occurs  because  our  previous  analysis  differed  from  theirs  only  by  our  allowing  total  pressure  to  vary  along  the  pore,  and 
by  consequences  of  that  variation.  Our  present  analysis  also  takes  into  account  that  the  fuel  concentration  is  nonzero  under  anode 
saturation  conditions,  unlike  our  previous  analysis  and  that  of  Jiang  and  Virkar.  The  result  as  seen  in  Table  1  is  that  our  present 
tortuosity  values  are  smaller  than  our  previous  values  by  roughly  10%.  We  consider  that  these  lower  values  are  more  accurate,  but 
it  must  be  remembered  that  the  previous  analyses  by  JV  and  ourselves  were  independent  of  a  model  for  V(i).  Our  confidence  on  the 
approximate  validity  of  this  10%  reduction  is  based  on  the  accuracy  with  which  our  V(i)  model  [8]  fits  the  V{i)  curves  of  JV  for  the 
H2/H2O  binary  system. 

Plots  of  gas  concentrations  vs.  position  in  the  anode  appear  in  Fig.  2  for  the  lowest  and  highest  H2  and  CO  fuel  concentrations 
in  the  anode  plenum.  The  total  concentration,  and  thus  total  pressure,  has  much  stronger  position  dependence  for  H2  than  for  CO 
fuel.  This  strong  dependence  of  total  pressure  on  position  in  the  anode  shows  that  the  commonly  made  assumption  of  uniform  total 
pressure  across  the  anode  is  a  serious  approximation.  The  quadratic  contribution  to  the  position  dependences  for  the  various  fuel 
and  exhaust  concentrations  is  barely  visible. 


6.  Calculation  for  ternary  gas  input 

The  ternary  systems  consisted  of  incoming  H2  fuel  diluted  with  He,  N2,  or  CO2.  They  are  ternary  because  H2O  (steam)  “exhaust” 
is  produced  at  the  anode-electrolyte  interface.  We  will  compare  the  JV  calculations  of  tortuosity  r  with  values  found  from  our 
analysis. 

In  the  preceding  analysis  of  the  binary  systems,  we  used  subscripts  e  and  f  for  exhaust  and  fuel  gases,  respectively.  For  the  ternary 
systems,  we  will  use  subscripts  d  for  the  He,  N2,  or  CO2  diluent  gases,  e  for  the  H2O  exhaust  gas,  and  f  for  the  H2  fuel  gas. 

We  note  that  for  H2  fuel  and  H2O  exhaust,  in  the  Stefan-Max  well  equations  the  molar  flow  rates  obey  Nd  =  0  and  Ve  =  —Nf  and 
we  denote  Nf  by  N.  Then  Eq.  (1)  becomes 

N  XdN  (Xf  +  Xe)N  _ 1  dQPZf) 

T>Kf  T>df  Def  RT  dx 

N  _  XdN  _  (Xf  +  Xe)N  =  1  d (PXe) 

Dkq  Dde  T>ef  RT  dx 

_XdN  +  Xd^  = _ 1  d  (PXd) 

D  df  £>de  RT  dx 

The  third  equation,  for  the  diluent  gas,  is  absent  for  binary  systems. 

First  we  convert  to  molecular  units  by  multiplying  through  by  Avogadro’s  number.  This  converts  N  to  /,  which  has  units  of 
molecules  cm-2  s-1,  and  converts  R  to  k  which  is  Boltzmann’s  constant.  Thus  URT  becomes  1/kT.  From  the  ideal  gas  law,  P  =  nkT , 
so  (l/kT)d(PXf)/dx  =  d(nXf)/dx  =  drif/dx,  etc.  These  are  the  concentration  gradients  in  cm-4.  Here  n  is  the  total  gas  concentration  and 
rif ,  etc.  are  the  partial  concentrations,  in  cm-3. 

We  can  then  write  Eq.  (25)  in  terms  of  molecular  units  as 


dnf 

J 

JX  d 

J(Xe  +  Xf ) 

dx 

Dk\ 

Ddf 

Dg  f 

d  n& 

J 

D  Ke 

JX  d 

Dde 

J(Xe  +  Xf) 

dx 

Def 

dnd 

JX  d 

JXd 

dx 

Ddf 

Dde 

The  quoted  D#  values  in  JV  Table  1,  which  we  call  D[j  1,  are  based  on  1  atmosphere  pressure,  which  we  label  P\.  Accordingly, 
using  the  ideal  gas  law  P  =  nkT ,  noting  that  P  is  in  the  denominator  of  Dy,  and  letting  n\  be  the  total  gas  concentration  at  T-  1073  K, 

1  _  P  _  n 

Dij  P\  Di j  1  n\Dij\ 


(27) 
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Finally,  we  note  that  X\  =  njn ,  etc.,  so  we  can  rewrite  Eq.  (26)  as 


drif 
dx 
d  ne 
dx 


J  J  \  n& 

£>Kf  ^1  |_Alfl 


J  J 

£>Ke  ^1 


Alel 


+ 


(ne  +  n/) 

Ajfi 
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T>efl 


d^d  _  ( Jnd\  f  1  1  \ 

dx  V  hi  /  V  Alel  Alfl/ 

Adding  these  three  equations  and  noting  that  rid  +  nQ  +  rif  =  n  yields  an  equation  for  n  and  its  solution, 


d  n  (l  1  \ 

—  =  •/  7. - 7^—  ,  n  =  n  i(l  +  c5x);  c5 

dx  \DKe  DKfJ 


In  solving  the  dn^/dx  equation  in  Eq.  (28)  for  n&,  we  get 


n<d  =  ^dp  exp  (c7x); 


<  0. 


(28) 


(29) 


(30) 


Jiang  and  Virkar  [1]  pointed  out  this  exponential  nature  of  the  diluent  gas  concentration  dependence  on  position,  but  omitted  using 
this  dependence  to  simplify  their  calculations.  This  non- vanishing  concentration  dependence  on  x  may  seem  surprising  because  the 
diluent  gas  has  no  net  flow  in  the  pores.  It  occurs  because  of  the  unlike  masses  of  the  fuel  and  exhaust  molecules  that  the  diluent 
molecules  collide  with.  The  >0  and  <0  relations  occur  in  Eqs.  (29)  and  (30),  respectively  because  for  the  systems  studied,  the  exhaust 
molecule  is  more  massive  than  the  fuel  molecule  and  the  inverse  Knudsen  and  binary  diffusion  coefficients  increase  with  molecular 
mass.  Accordingly,  Eq.  (29)  tells  us  that  total  gas  concentration  increases  linearly  with  x,  whereas  Eq.  (30)  tells  us  that  diluent  gas 
concentration  decreases  exponentially  with  x. 

Because  n\  and  nap  are  known  as  explained  below,  we  can  solve  for  nt  in  terms  of  rif 


—  n  —  rid  —  nf  =  n\(\  +  c5x)  —  /?dp  exp ( cjx )  —  rif. 


(31) 


Then  the  drif/dx  equation  in  Eq.  (28)  becomes 


drif 

dx 
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exp  (c7x). 


The  solution  to  this  equation  can  be  written  as 


nf  =  nfp  +  ndpg[l  exp  (c7x)] 


x  — 


J 

T^efl 


(x  +  -c5xz); 
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(32) 


(33) 


Eq.  (33)  is  especially  useful  under  the  condition  of  anode  saturation.  Under  that  condition  we  previously  [5]  set  rif  equal  to  zero  at 
the  electrolyte  interface,  where  x  =  L  =  zw ,  but  as  explained  in  the  binary  system  analysis,  we  now  set  rif  =  ftfas,  its  nonzero  value 
for  the  current  density  value  ias  for  which  the  output  voltage  V  drops  rapidly  through  zero  on  its  way  to  —00. 

We  also  need  expressions  for  ridv  and  rifv  under  the  anode  saturation  condition.  To  find  these,  we  note  that  under  this  condition 
the  net  flow  rates  in  molec  s-1  into  the  anode  plenum  are 


jd  =  jiPd, 


jf  =  71  (1  Pd) 


iasS 

q 


(34) 


The  total  of  these  flow  rates  into  the  plenum  is  simply  j\ .  Here,  ias  is  the  anode  saturation  current  density  as  given  in  JV  Table  2, 
S  is  the  electrode  area  of  1.1  cm2,  q  is  the  charge  of  3.2  x  10“ 19  C  per  gas  molecule  created  or  destroyed  at  the  anode/electrolyte 
interface,  and pd  is  the  fractional  partial  pressure  of  the  diluent  gas,  given  in  JV  Table  2  as  Pn2  > or  Pco2  •  Note  that  the  superscript 

o  used  in  Eq.  (26)  of  JV  for  metered  incoming  gas  is  missing  in  JV  Table  2.  One  finds  the  total  incoming  molecular  flow  rate  j\  from 
the  ideal  gas  law,  Pi  V[n  =  NmkTm,  using 


dMn  =  Pi  dUin 

dt  kTm  dt 


(35) 


using  dVin/d t  as  140  cm3  min  1  as  stated  by  JV  for  all  experiments.  In  Eq.  (35)  we  use  Tm  =  273  K  to  best  fit  the  results  of  JV  using 
their  calculational  method,  so  j\  =  6.286  x  1019  molec  s-1. 
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As  discussed  for  the  binary  systems,  the  net  flow  rates  into  the  plenum  equal  the  flow  rates  out  of  the  plenum  exit  orifice,  and 
these  exit  orifice  flow  rates  are  proportional  to  the  gas  concentrations  in  the  plenum,  so 


J  d 

ft  dp  =  — ftl  =  Pdftl, 
J 1 


J  f 

ftfp  =  —ftl 
Jl 


1  ks^\  ( \  ks\  .  Qj\ 

i  -  Pd —  Ui  =  1  -  Pd  -  —  J  fti;  n  =  — . 

qji  J  \  n )  s 


Finally,  we  find  J  under  the  anode  saturation  condition,  as  discussed  for  the  binary  systems 


Jas  — 


<Pq 


(36) 


(37) 


We  can  now  write  Eq.  (33)  under  the  anode  saturation  condition  and  evaluated  at  v  =  ra>.  From  Eqs.  (29),  (30),  and  (37),  we  see 
that  C5  and  C7  depend  linearly  on  t.  Accordingly,  we  must  use  the  r-independent  factors  C6  =  c$w/z  and  c8  =  C'jw/r  to  display  the 
true  r  dependence. 


ftfasC*  =  L  =  rw)=  1  -  Pd  ~  —  ftl  +  Pdn\g[l  -  exp  (c8r2)]  - 
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(38) 


This  is  a  transcendental  equation  for  r2  that  can  be  solved  numerically.  Written  out  explicitly  by  substituting  for  C5,  C7,  and  g 
from  Eqs.  (23),  (24),  and  (27),  respectively,  it  appears  as 
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In  Eq.  (39),  we  use  ftfas  =  fti4s/ftoi  as  we  did  for  the  binary  system. 


2  \(j)qn\  J  V^Ke  T>Kf/ 

0. 


(39) 


7.  Results  for  ternary  gas  input 

The  relations  in  the  equations  below  are  useful  in  solving  Eq.  (39)  numerically  for  the  tortuosity  t  values  presented  in  Table  2 
and  Fig.  3.  The  Z)de  1  values  below  are  calculated  from  Eq.  (2.22)  and  Tables  2.2  and  2.3  of  Mass  Transfer  [7].  The  constants  that 
remain  the  same  from  one  run  to  the  next  are 


Fig.  3.  Plots  of  gas  concentrations  vs.  position  in  anode  for  some  ternary  gas  inputs  under  anode  limiting  current  conditions. 
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£>Ke  =  3-767,  £>Kf  =11.3,  £>efi  =  7.704  cm2  s_1, 
Z)dei  =  7.462,  £>dfi  =  13.29  for  d  =  He, 

^dei  =  2.21,  Z)dfi  =  6.303  for  d  =  N2, 

Ddei  =  1.73,  Dd fi  =  5.56  for  d  =  C02, 
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=  0.05877  s  cm-2  for  d  =  He, 
=  0.2938  for  d  =  N2, 

=  0.3982  for  d  =  C02, 

0.0546  s  cm-2  for  d  =  He, 

=  -0.0289  for  d  =  N2, 

=  -0.0501  for  d  =  C02, 


£>Kf  Ajfl 


=  0.2183  s  cm 
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=  0.02297  s2  cm-4, 


(40) 


Pi  =  1.015  x  105N  m 


-2 


k  =  1.38  x  1(T23JK, 


T  =  1073  K, 


w  =  0.11  cm, 


5  =  1.1cm2, 
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n\  =  —  =  6.855  x  1018cm  3, 
kT 

=  0.54,  q  =  3.2  x  1(T19  C,  ji  =  6.286  x  1019  s"1, 


aoi  =  48.4  A  cm  2,  i\  =  —  =  18.29  Acm  2, 
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5.634  x  1016, 


2(j>qn\  \  DKe  Acf/ 


=  0.008217  cm2sC_1, 


4>qni 

\  T^del 

Ddi\ 

w  ( 

(  1 

1 

<pqn\ 

\  T^del 

T^dfl 

w  ( 

(  1 

1 

(jiqni 

\  T^del 

T^dfl 

=  0.005458  cm2s  C"1  for  d  =  He, 


=  0.02728  cm2s  C_1  for  d  =  N2, 


=  0.03698  cm2s  C_1  for  d  =  C02. 


Though  Eq.  (38)  is  easy  to  solve  exactly  for  t  by  numerical  methods,  a  sufficiently  accurate  analytical  expression  for  r  can  be 
obtained  by  approximating  the  exponential  term  by  a  quadratic  expression, 

exp  (c8t2)  =  1  +  c8t2  +  2(c8)2t4.  (41) 

Then  Eq.  (38),  or  Eq.  (39),  becomes  a  quadratic  equation  in  r2  and  can  be  solved  for  r  in  the  usual  manner.  This  quadratic  equation, 
its  solution,  and  its  Ar,  Br  and  C  coefficients  are 
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If  we  ignore  the  x2  and  thus  the  r4  terms  in  the  foregoing  equations  by  setting  A '  =  0,  we  find 


7  \  1/2 


B' 


^linear  — 


(45) 
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Regarding  the  zb  sign  choice  in  Eq.  (42),  we  note  that  B'  is  negative  and  C'  is  positive.  If  A!  is  positive,  we  must  choose  the  + 
sign  because  it  gives  the  correct  limit  given  in  Eq.  (45)  as  A!  goes  to  zero.  If  A!  is  negative,  we  must  choose  the  —  sign  to  avoid 
taking  the  square  root  of  a  negative  number  in  Eq.  (42).  In  Eq.  (43),  all  expressions  in  round  parentheses  are  positive,  except  that 
(1/Z)efi  —  l/Z^dfi )  is  negative  for  N2  and  CO2  diluent  gases.  Thus  A!  is  negative  for  He,  but  for  N2  and  CO2  diluent  gases  the  sign  of 
A!  depends  on  which  term  dominates  in  Eq.  (43). 

Results  for  the  ternary  gas  systems  are  displayed  in  Table  2  below.  There  is  a  systematic  decrease  in  the  calculated  tortuosity  as 
the  diluent  concentration  decreases,  except  for  our  calculations  for  He.  Also,  there  is  a  systematic  increase  in  calculated  tortuosity 
as  diluent  molecular  mass  increases.  In  theory,  the  tortuosity  should  be  independent  of  such  experimental  conditions.  Although  the 
reasons  for  these  trends  are  not  known,  we  do  obtain  the  expected  result  that  the  t  values  based  on  the  three  diluents  approach  each 
other  as  diluent  concentration  decreases. 

The  linear  approximation  to  r  is  accurate  within  a  few  percent.  The  r  values  calculated  by  JV  (the  square  roots  of  the  “tjv” 
values  reported  in  their  paper,  which  are  really  “tortuosity  factors”)  are  larger  than  ours  for  the  light  He  diluent  gas,  but  comparable 
to  ours  for  the  heavy  N2  and  CO2  diluent  gases.  The  variation  in  our  “rexact”  values,  and  the  differences  between  our  “Texact”  values 
and  the  corresponding  (7jy)1/2  values,  is  greater  than  for  the  binary  systems  presented  in  Table  1.  Still,  all  of  our  reXact  values, 
centered  near  2.3  and  ranging  from  about  1.7  to  3.3,  are  physically  reasonable  and  provide  confidence  that  JV’s  experimental 
method  is  valid  for  finding  tortuosity. 

Plots  of  gas  concentrations  across  the  anode  for  four  ternary  systems  under  anode  limiting  current  conditions  appear  in  Fig.  3.  All 
show  significant  increases  in  total  gas  concentration  n  (and  thus  in  total  pressure  P )  across  the  anode.  All  show  decrease  in  diluent 
concentration  going  across  the  anode,  because  collisions  with  the  heavier  outflowing  H2O  molecules  have  more  effect  than  those  of 
the  lighter  inflowing  H2  molecules,  even  though  the  H2O  and  H2  molecular  flow  rates  are  the  same.  The  exponential  curvature  of 
the  diluent  concentrations  is  small  but  noticeable,  whereas  the  quadratic  and  exponential  corrections  to  the  linear  v  dependences  of 
fuel  and  exhaust  concentrations  appear  insignificant. 

8.  Conclusions 

Based  on  all  the  above  results,  we  consider  that  the  tortuosity  for  the  Jiang  and  Virkar  apparatus  [1]  analyzed  herein  is  2.3  zb  0.6. 
This  is  fairly  consistent  with  the  range  of  2.5-3  found  by  Williford  et  al.  [3]  for  modem  SOFC  materials,  and  is  also  within  the  1-3 
range  that  seems  most  physically  reasonable,  1  being  the  limit  for  an  extremely  open  and  probably  impractical  electrode  design.  The 
apparent  disagreement  between  our  values  and  the  tortuosity  values  reported  by  Jiang  and  Virkar  lies  in  their  reporting  as  tortuosity 
what  is  often  called  the  tortuosity  factor,  which  is  the  square  of  the  tortuosity  [2]  so  we  can  say  that  their  actual  values  and  ours  and 
those  of  Williford  et  al.  are  in  general  agreement  and  are  physically  reasonable. 

The  more  exact  method  presented  here  for  calculating  tortuosity,  or  alternatively  for  calculating  concentration  drops  across  the 
electrode  if  tortuosity  is  known,  is  simple  enough  to  justify  its  use.  Knowing  concentrations  at  the  electrode-electrolyte  interfaces  is 
essential  for  calculating  the  concentration  polarization  term  in  the  voltage  vs\  current  density  expression.  Conversely,  our  expressions 
for  activation  and  concentration  polarizations  based  on  opposing  current  flows  at  each  electrode  interface  with  the  electrolyte 
presented  herein  and  derived  elsewhere  [8]  allow  us  to  determine  the  tortuosity  more  accurately.  Such  voltage  calculations,  together 
with  calculations  of  pressure  drops  across  electrodes  and  of  tortuosity  as  presented  above,  provide  tools  for  predicting  performance 
of  new  SOFC  designs  and  for  analyzing  performance  of  existing  designs,  both  in  the  fuel  cell  and  electrolysis  operation  modes. 
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